Lattice Ising model in a field: Eg scattering theory 

V. V. Bazhanov*f, B. Nienhuis* and S. O. Warnaar^ 



(N 



CO 
On 
G\ . 

Abstract 

° . 

Zamolodchikov found an integrable field theory related to the Lie algebra 
Q ' Eg, which describes the scaling limit of the Ising model in a magnetic field. 

He conjectured that there also exist solvable lattice models based on Es in 
the universality class of the Ising model in a field. The dilute A3 model is a 
solvable lattice model with a critical point in the Ising universality class. The 
parameter by which the model can be taken away from the critical point acts 
like a magnetic field by breaking the W.2 symmetry between the states. The 
expected direct relation of the model with Es has not been found hitherto. In 
this letter we study the thermodynamics of the dilute A3 model and show that 
in the scaling limit it exhibits an appropriate Es structure, which naturally 
Q\ [ leads to the Es scattering theory for massive excitations over the ground state. 

Ph! 1 Introduction 

r-j 

Since the work || by A.B. Zamolodchikov it is known that certain perturbations 
of conformal field theories (CFT's) lead to completely integrable models of massive 
quantum field theory (QFT). The existence of non-trivial higher integrals of motion 
and other dynamical symmetries [§, || §, [|, || in such a QFT allows to compute 
the spectrum of the particles and their 5"-matrix explicitly. At the same time, these 
QFT models can be obtained as the scaling limit of appropriate non-critical solvable 
lattice models in statistical mechanics (see [0] for an introduction and references on 
solvable lattice models). In the latter approach the spectrum and the S"-matrices can 
be calculated from the Bethe Ansatz equations for the corresponding lattice model 
p. |9|, [10]. The natural problem arising in this connection is to find lattice counterparts 



for all known integrable perturbed CFT's and vice versa. A description of known 
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results of such correspondence lies outside the scope of this letter and we refer the 
interested reader to [1-10] and references therein. Here we consider one particularly 
important example of this correspondence associated with the Ising model at its 
critical temperature in a magnetic field, hereafter referred to as the magnetic Ising 
model. 



A.B. Zamolodchikov has shown [11] that the c = 1/2 CFT (corresponding to 
the critical Ising model) perturbed with the spin operator 12 = 02,2 of dimension 
(1/16,1/16) describes an exactly integrable QFT containing eight massive particles 
with a reflectionless factorised S'-matrix. Up to normalisation the masses of these par- 
ticles coincide with the components Sj of the Perron- Frobenius vector of the Cartan 
matrix of the Lie algebra E 8 : 

- = f. (LI) 

rrij Sj 

The element of the S'-matrix describing the scattering of the lightest particles, with 



mass mi, reads 11 



th (f- if) th(f- if) th(f- if;)' (L2) 



with (3 the rapidity. The other elements are uniquely determined by the bootstrap 
program p"T|j . 

The aim of this letter is to show that the above QFT describes the scaling limit 
of the dilute A3 model of Warnaar, Nienhuis and Seaton |12|, [13[ in the appropriate 



regime. It should be noted that there were some earlier, rather strong indications 
supporting the above correspondence. All these parts remarkably fit together with 
our results, completing a sequence of arguments which can be summarised as follows: 

(i) The dilute A3 model is an interaction-round-a-face model on the square lattice 

with spins taking three values (detailed definitions are given in equations ( [2.1|) - 
( |2.4| )). Admissible values of the adjacent spins are determined by the incidence 
matrix (|2.1|) , which has largest eigenvalue equal to 1 + 

(ii) The model has two physically distinct regimes of relevance to our discussion, 

here denoted as i) and ii), depending on the region of the spectral parameter 
or, equivalently, of a sign of the Hamiltonian of the associated one-dimensional 



chain. (These are the regimes 2 + and 3 + of ref [13], respectively). The central 
charges and the conformal dimensions of the leading perturbation computed 
from exact expressions for the free energy and the local state probabilities of 



the dilute A 3 model for these two regimes read [|14], [13[ 



i) c=l/2, A = 1/16; ii) c = 6/5, A = 15/16. (1.3) 

(iii) In ref |J pl| Bazhanov and Reshetikhin proposed thermodynamic Bethe Ansatz 
equations (TBAE) related to the A-D-E Lie algebras, corresponding to non- 
critical models in statistical mechanics. Using standard thermodynamics cal- 
culations and the high level Bethe Ansatz (see M and references therein) they 



computed: the central charges of the corresponding scaling field theories, di- 
mensions of the leading perturbations, the spectra and scattering amplitudes 
of the massive excitations, expressing them through fused Boltzmann weights. 
In particular, in the case relevant to our discussion (Q—E s , g = 30, p = £ = 1, 
in the notation of H) the exponents they found Q precisely match ( |1.3| ) in both 



regimes. Furthermore, the TBAE allowing the calculation of the largest eigen- 
value of the incidence matrix of the underlying lattice model, gave in this case 
precisely the value 1 + v2 |16| . 

(iv) Finally, the spectrum and S'-matrix of the scaling field theory in regime i) found 
in H from the high level Bethe Ansatz for Eg coincide with those of Zamolod- 
chikov's magnetic Ising model. 

All the above arguments strongly suggest that the TBAE based on the Lie algebra 



E 8 as proposed in || [L5| , are those of the the dilute A 3 model. 

In this paper we present the Bethe Ansatz equations (BAE) for the non-critical, 
dilute A^ model. As these equations, at criticality, are very similar to those of the 



Izergin-Korepin model [T7|, |TB|, it is not surprising that, when specialised to L = 3, 
they do not display any explicit structure related to the root system of Eg. It turns 
out however that this structure reveals itself in a quite complicated string structure 
of the solutions to the BAE. Motivated by an extensive numerical investigation of the 
BAE we formulate an exact conjecture for the thermodynamically significant strings. 
This leads to TBAE, which, rewritten in a new string basis precisely yield the Eg 
based TBAE of ref || discussed under (iii). As a result of (iv) this finalises the 
correspondence between the dilute A3 model and the magnetic Ising model. 



2 The dilute A models 

The dilute Al model, belonging to the more general class of dilute A-D-E models, is 
an exactly solvable, restricted solid-on-solid model defined on the square lattice. Each 
site of the lattice can take one of L possible (height) values, subject to the restriction 
that neighbouring sites of the lattice either have the same height, or differ by ±1. This 
adjacency condition can be conveniently expressed by a so-called incidence matrix M: 

M afi = 5 0j6 _i + 5 a>b + 5 a>b+ i a,b G {1, . . . ,L}, (2.1) 

where we note that M relates to the Cartan matrix C Al of the Lie algebra A^ by 
M = 31 — C Al , with / the identity matrix. The eigenvalues of the incidence matrix 
are found to be 

A j = l + 2cos( z ^ T ) j = l,...,L. (2.2) 

1 Note that the equations (5-1) and (5-4) in || have been misprinted. Correcting (5-1) to c = 
c^{l) + <P(r — I — g) — <P (r — g) + rank Q yields the following result for the central charge in (5-4): 
c = 2 rank£/(<7 + 2). Also the phrases "minimal unitary" , just before, and "by the operator <^(i,3)" 
just after (5-4) should be deleted. 



For the case of interest here, L = 3, we thus find the largest eigenvalue to be 1 + \/2, 
in accordance with the prediction for the E 8 TBAE as mentioned in (iii) of the 
introduction. 

Using standard definitions of &i(u, q) -functions, suppressing the dependence on 
the nome q = e _T , r > 0, the Boltzmann weights of the allowed height configurations 
of an elementary face of the lattice are 
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The variable A and the range of the spectral parameter u in the above weights are 
given by0 

^ _ 7r L + 2 j < u < 3A regime i) , 

4L + 1 1 3A — 7r < u < regime ii). ^ ' ' 



2 In (T3] two more regimes were defined, which are omitted being of no relevance here. 



3 Bethe Ansatz 



The transfer matrix of the dilute A models is defined in the usual way as 



T, 



{b} 
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N 

l[w[ 

i=i 



Uj+l 

a j+ i 



where {a} is an admissible path of heights and cin+i 



(3.1) 

a-i, bN+i = b\. The number 



of admissible paths is given by Trace(A JV ) = £^ =1 Af • 



We define two positive integers p and r and a renormalised spectral parameter 9 



by 
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V 



e 



7T 
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r 2r 

where p and r are coprime. Furthermore we introduce the modified ^-function 



h{6) 



1 -i/4 Q H \ 



(3.2) 



(3.3) 



With these definitions the eigenvalues of the transfer matrix are found to be 
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(3.4) 



_J h(6) h(2p - g) \ N " h(6 + i 9j - Sp) 
+ U \ h{Ap) h{6p) J j-} 1 h{6+ iej-Ap) 

where the numbers {8j} are given by the following set of BAE: 



( h(2p+ 16,) ^ " hjiOj- i6 k +Ap) hjxdjj- iO k -2p) 
\h{2p-ie 3 )) I} 1 h(i9 j -i6 k -4p)h(i0 j -i9 k + 2p) 



J 



N (3.5) 



and uj = exp(i7r£/(L + 1)), i = 1,...,L. Note that if we set the nome q of the 
function h equal to zero and choose uj to be unity, the above eigenvalue expression 
and BAE reduce to those of the Izergin-Korepin model [18] in the sector in which the 
number of roots is equal to the system size N. A proof of equations (|3.4|) and (|3.5| ) in 
the critical limit, based on an extension of the mapping of RSOS models onto vertex 
models found in fl9 ], will be presented elsewhere [20|. In the general elliptic case the 
equations can be proven using the functional equations for the dilute A models as 



obtained in ref [21 



4 Thermodynamic Bethe Ansatz 



The remaining part of this letter is devoted to the solution of the BAE in the thermo- 
dynamic limit for the dilute A 3 model, corresponding to the case p = 5 and r = 16. 

Let n® be a positive integer, A® an n^-dimensional vector with integer coef- 
ficients A® and e® = 0, 1. We then define a string of type t as a set of complex 
numbers {o^} with 

a® = a® + i (A® + e®r) k = 1, . . . , n®, (4.1) 

where the real number a® specifies the centre of the string. Based on an extensive 
numerical study [g0| we find that for N — > oo each solution {6j} of fl3.5|) consists of 
a collection of strings, where we have only nine thermodynamically significant string 
types (in the sense that /N is finite for N — > oo) 

8 

AT = ^ n ®N® +o{N), (4.2) 

4=0 

with N® the total number of strings of type t. These nine allowed string types 
are listed in Table 1. We expect the above statement to be exact and claim it as a 
conjecture. 
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Table 1: The nine thermodynamically significant string types 

In the thermodynamic limit the centers of the strings form continuous distribu- 
tions and the BAE (|3.5|) lead to integral equations for the densities of strings pt{pi) 
and "holes" p t {a) 

8 

b t (a) = -(-) 5t -°p t (a) + Y<Bt, s *Ps(a) t = 0,...,8, (4.3) 
where a * b denotes the convolution of the functions a and b 

/rr /n 
a(a- (3)b{(3)d(3. (4.4) 
-tt/tt 



The functions b t and B t s in ( |4.3|) , which are 2rr/7r-periodic, read 

b t (a) = ^^(a+iAjj? 
k=i 

B t , s (a) = -(-)**M(«) + EE K (a + i A? + i A t w 



fc=l f=l 
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^ 2p (a+ iA^ + iA 



with ^ defined as 
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The function tp^ has the following Fourier transform (FT) 
sinh(r — k)x 



where 



F(x) 
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< k < 2r, 
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e" iaa; F(a) da 



2vr 



£ e iQ ^F(x n ), 



(4.8) 



with 5 = 7r 2 /(rr), x„, = 5n. 

As usual, we define the (local) Hamiltonian 7i of the associated one-dimensional 
integrable spin chain as the logarithmic derivative of the transfer matrix at 9 = 0. 
Then, after appropriate normalisation and shift, the spectrum of TC, in the limit 
N — > oo reads 



Tr = e E / TrA M«) /°*(«) da + ° ( e ~ MiV ) ^ > °- 

IV J —rr/n 



(4.9) 



where e = — 1 for regime and e = 1 for regime ii) in 
The densities p t are normalised such that 



tt/tt 

-TV 1 7T 



^(a) da = iV w /iV. 



Therefore, from equation ( |4.2j ), we have 

^ rrr/n 

E / n ■ Pt( a ) da = 1. 

This relation together with equation (|4.3|) for £ = implies 
p (a) = 0. 



(4.10) 



(4.11) 



(4.12) 



Hence we conclude that the strings of type have no holes in any state, and we 
eliminate po(a) from (|4.3| ). After a tedious calculation we find that the resulting 
integral equations can naturally be described in terms of the E 8 root system as follows. 

Let Cfg t, s — 1, . . . , 8 be the elements of the Cartan matrix for Eg, where we use 
the following enumeration of the nodes of the corresponding Dynkin diagram: 

.8 



1 2 3 4 5 6 7 

Furthermore, define the functions Kfj, Afl, af| and s by their FT 

= 5 t , s + s(x)(cf; s -25 t , s ) 
ig(x) = [k Es (x)]^ 

ag(x) = (4.13) 

six) = 

2 cosh x 

With these definitions, and after eliminating p , the integral equations ( [4.3| ) and the 
energy expression Q4.9| ) take the form 



ag(a) = p t (a) + J2 A ?, 8 s*Ps(a) t = l,...,8, 

s=l 

E 8 rrr/n 

— = -e^2 / af 8 t (a) p t (a) da + const. (4-14) 



t=l 



We can now use ( |4.14[ ) to study the scaling limit of the model. In fact, all relevant 
calculations have already been carried out in ref [[J and we only need to refer to the 
appropriate results therein. To make the correspondence with ref somewhat more 
transparent, let us give the expression for the equilibrium free energy F(T) of the one- 
dimensional spin chain at finite temperature T, as it follows from ( [4.14j ) via standard 



TBA calculations [22 



F(T) 8 f TT / 7T 
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/TV I TV , . 

a\\{a) Tlog (1 + e-^'W) da + const, (4.15) 



where j3 = 1/T is the inverse temperature. The functions e t = T\og(p t / p t ) are the 
solutions of the integral equation 

8 

eSi >t s(a) =Tlog(l + e- /3et(Q) ) - * Tlog (l + e^ s(a) ) (a). (4.16) 

s=l 

The above two equations are equivalent to (3-20) and (3-21) of ref 0], respectively, 
with their ^=E§, r = 32, g = 30, p — £ — 1, their nome q replaced by q 1 ^ 2 and with 
their e" negated. This last difference reflects the fact that our TBA equations are dual 



to those of ref |J in the sence that the densities of strings and holes are interchanged. 
From fl4l5p and (gig) it follows that for T = 

i) e = -1 et(a) = a§ 

(4.17) 

ii) e = +1 e t (a) = — 5 t is(a). 



The functions |et(a;) | are the energies of the excitations over the ground state. 

For e = — 1 the ground state is formed by type strings. As was remarked 
after equation ([4.12 ), these strings have no holes for any state. Therefore the Dirac 



sea is "frozen", and the excitations correspond to the remaining eight string types. 
The phenomenon of "freezing" of the Dirac sea which can be interpreted as the 



confinement of "holes" has been first observed in the TBAE calculations of ref IF23 



for the RSOS models of Andrews, Baxter and Forrester [24 



For e = 1 the Dirac sea is formed by the type 1 strings, and the only excitations 
correspond to holes in the Dirac sea. These excitations are of the kink type. 

Now we consider the scaling limit. We introduce a dimensional spacing parameter 
d for our chain and take the limit iV — > oo, d — *■ 0, keeping the (dimensional) 
length of the chain L = Nd to be macroscopically bigger than the correlation length: 
L » R c = q^^d, where £ is the index of the correlation length. In the scaling limit 
we thus have d ~ q^, N >> q~^, q — > 0, and we obtain the massive relativistic 
spectrum of excitations. To find this, one has to compute the energy dispersion law 
for the physical excitations in the q —>■ limit keeping the rapidities a of the order 
of «o — tt/tt, where the functions |e 4 (a)| have their minima. Taking into account 
the correspondence in notation discussed after (|4.16|) , one gets from (4-1) and (4-2) 
of ref 

i) e t(fP + oc Q ) = m t cosh (3 + otf) 
m t = const S t q*, £ = ^ 

ei(f/3 + a ) = mcosh/? + o(g ? ) (4.18) 
m = const q^, £ = 8, 



where (3 here denotes the rapidity variable and S t was defined just before equation 
(©■ 

Using the scaling relation £ = (2 — 2A) 1 it is seen that the values of £ in ( |4.18[ ) 
lead exactly to the dimensions of the leading perturbations as given in (|1.3|) . 

The values of the central charges of the corresponding (ultraviolet) conformal field 



theories listed in ( |1.3| ) have also been previously calculated. For regime i) in 0. [25]. |26 
and for regime ii) in ||. 

Finally, the S'-matrix for regime i), where all the string excitations for t = 1 ... 8 
correspond to distinct particles, can be found straightforwardly from equation ( 4.14Q . 
The result is 

^^pJifAgM ^ fcl (4.19) 



which coincides with Zamolodichkov's E§ S'-matrix 



For regime ii) the kink-kink S'-matrix is of the RSOS type related to the E 7 Lie 
algebra ||, and will be discussed elsewhere [pOR . 



5 Summary and Conclusion 

In this paper we have established the final link between Zamolodchikov's E 8 S-matrix 



of the critical Ising model in a field [|11J] and its underlying lattice model. By making 
a conjecture for the possible string solutions of its Bethe Ansatz equations, we have 
derived a system of thermodynamic BAE for the dilute A 3 lattice model of Warnaar 



et al. fl2 |. After a suitable transformation we have recast these TBAE in terms of 
the root system of the Lie algebra Eg. These Eg TBAE are found to be precisely 
those conjectured earlier by Bazhanov and Reshetikhin ||, and using their results, 
the correspondence between the dilute A 3 model and the Eg S'-matrix is made. 

To conclude we mention that two more remarkable integrable 0^2 perturbations 
of CFT's are known, notably those related to S-matrices with hidden E 7 (c=7/10) 
and E 6 (c=6/7) structure f27fl . Like the E 8 case, the underlying lattice models of 
these integrable QFT's correspond to models in the dilute A hierarchy. The working 
for these two extra cases, corresponding to dilute A 4 and A 6 , respectively, as well 
as some additional results for the dilute A 3 model will be the subject of a future 
publication 
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